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In the spin-excitation-mediated pairing mechanism for superconductivity, the geometric frustra¬ 
tion effects not only the spin configuration but also the superconducting-state properties. Within 
the framework of the kinetic-energy-driven superconducting mechanism, the doping and tempera¬ 
ture dependence of the Meissner effect in triangular-lattice superconductors is investigated. It is 
shown that the magnetic-field-penetration depth exhibits an exponential temperature dependence 
due to the absence of the d-wave gap nodes at the Fermi surface. However, in analogy to the dome¬ 
like shape of the doping dependence of the superconducting transition temperature, the superfluid 
density increases with increasing doping in the lower doped regime, and reaches a maximum around 
the critical doping, then decreases in the higher doped regime. 

PACS numbers: 74.25.N-, 74.20.Mn, 74.20.Rp, 74.72.Ek 


I. INTRODUCTION 


Investigation of oxide compounds has uncovered many 
unusual properties characterized by the strong electron 
correlation, which include unconventional superconduc¬ 
tivity and anomalous properties in the normal-state^. 
Superconductivity in cuprate superconductors results 
from some special microscopic conditions^ - — : (a) the one- 
half spin Cu ions situated in a square-planar arrange¬ 
ment and bridged by oxygen ions; (b) the weak cou¬ 
pling between neighboring layers; and (c) the charge- 
carrier doping in such a way that the Fermi level lies 
near the middle of the Cu-0 a* bond. One common 
feature of cuprate superconductors is the square-planar 
Cu arrangement^ - —. However, some oxide materials with 
a two-dimensional spin arrangements on non-square lat¬ 
tices have been synthetized^^. In particular, it has been 
reported^ - — that there is a class of cobaltate supercon¬ 
ductors Na^CoC^ • 2 /H 2 O, which have a lamellar struc¬ 
ture consisting of the two-dimensional C 0 O 2 layers sep¬ 
arated by a thick insulating layer of Na + ions and H 2 O 
molecules, where the one-half spin Co 4+ ions sites sit 
not on a square-planar, but on a triangular-planar lat¬ 
tice, therefore allowing a test of the effect of the geo¬ 
metric frustration on superconductivity^ - —. Moreover, 
Na^CoC^ • yH 2 0 is viewed as an electron-doped Mott in¬ 
sulator, where superconductivity appears with electron 
doping^ - —. Furthermore, it has been found that the an¬ 
tiferromagnetic (AF) spin correlation- 41-15 is present in 
Na^-CoCVy^O, although being much weaker than those 
in square-lattice cuprate superconductors. In this case, 
a question is whether the unusual features observed on 
square-lattice superconductors exist also in triangular- 
lattice superconductors or not? The finding of supercon¬ 
ductivity in triangular-lattice cobaltate superconductors 


has raised the hope that it may help solve the unusual 
physics in square-lattice cuprate superconductors. On 
the other hand, the doped Mott insulator on a triangu¬ 
lar lattice is also of interests in its own right with many 
unanswered fascinating questions^d^, where the geomet¬ 
ric frustration was expected to destroy the AF long-range 
order (AFLRO) and leads to a quantum spin-liquid state. 

Superconductivity is characterized by exactly zero 
electrical resistance and expulsion of magnetic fields oc¬ 
curring in superconductors when cooled below T c . The 
later phenomenon is so-called Meissner effect^, i.e., a 
superconductor is placed in an external magnetic field B 
smaller than the upper critical field B c , the magnetic 
field B penetrates only to a penetration depth A and 
is excluded from the main body of the system. This 
magnetic-field-penetration depth is a fundamental pa¬ 
rameter of superconductors, and provides a rather direct 
measurement of the superfluid density p s (p s = A -2 )i£, 
which is proportional to the squared amplitude of the 
macroscopic wave function. In particular, the varia¬ 
tion of the magnetic-field-penetration depth as a function 
of doping and temperature gives the information about 
the nature of quasiparticle excitations and their dy¬ 
namics. Moreover, the magnetic-field-penetration depth 
can be also used as a probe of the pairing symmetry, 
since it can distinguish between a fully gapped and a 
nodal quasiparticle excitation spectrum^ -19 . The for¬ 
mer results in the thermally activated (exponential) tem¬ 
perature dependence of the magnetic-field-penetration 
depth, whereas the latter one implies a power law be¬ 
havior. For square-lattice superconductors, the Meiss¬ 
ner effect has been studied experimenta ll y ^ S' ■ - 1 as well 
as theoreticall y ^. 2 1 In particular, the electromag¬ 
netic response in square-lattice superconductors has been 
discussed^ based on the kinetic-energy-driven supercon¬ 
ducting (SC) mechanis m 24 ! 25 , and the obtained results of 
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the doping and temperature dependence of the magnetic- 
field-penetration depth and superfluid density are well 
consistent with the experimental data observed on 
square-lattice superconductor s 20 ’ 21 . In triangular-lattice 
cobaltate superconductors, on the other hand, although 
the Meissner effect has been investigated by virtue of sys¬ 
tematic studies using the muon-spin-rotation measure¬ 
ment techniqu e) 26 ’ 27 , the electromagnetic response has 
not been clarified starting from a microscopic SC theory, 
and no explicit calculations of the evolution of the super- 
fluid density with doping and temperature has been made 
so far. In this paper, we try to study this issue within the 
framework of the kinetic-energy-driven SC mechanism. 
We show that the magnetic-field-penetration depth of 
triangular-lattice superconductors exhibits an exponen¬ 
tial temperature dependence due to the absence of the 
d-wave gap nodes. However, in analogy to the case of 
square-lattice superconductors, the superfluid density in 
triangular-lattice superconductors also has a dome-like 
shape of the doping dependence. 

The rest of this paper is organized as follows. The basic 
formalism is presented in section El where we generalize 
the response kernel function obtained within the frame¬ 
work of the kinetic-energy-driven SC mechanism from the 
case in the previous square-lattice superconductors^ to 
the present case for triangular-lattice superconductors, 
and then employ this response kernel function to obtain 
explicitly the doping dependence of the Meissner effect 
in triangular-lattice superconductors for all the temper¬ 
ature T < T c . Based on this theoretical framework of 
the electromagnetic response, we then discuss the basic 
behavior of triangular-lattice superconductors in a weak 
electromagnetic field in section IIIII Finally, we give a 
summary in section m 

II. THEORETICAL FRAMEWORK 

In triangular-lattice cobaltate superconductors, the 
common feature is the presence of the C 0 O 2 planed—, 
and then it is thus believed that the nonconventional SC 
mechanism in triangular-lattice cobaltate superconduc¬ 
tors and the related anomalous properties in the normal- 
state are dominated by this plane. In this case, many 
authors have argued that the essential physics of the 
C 0 O 2 plane is contained in the t-J model on a trian¬ 
gular lattice^. To study the electromagnetic response 
in triangular-lattice cobaltate superconductors, the t-J 
model can be generalized by including the exponential 
Peierls factors as, 

H = i^ e_i(e/B)A( * HpC ^ Ci +^ pt 

lfj<7 

- i‘T. pc l c ‘- p, + J Y. s l ■ Sl+f/, (1) 

la If) 

where the electron hopping integral t < 0, the summation 
is over all sites l , and for each l, over its nearest-neighbor 


7 ), Cj a and Cia are operators that respectively create and 
annihilate electrons with spin cr, S( = (Sf, Sf , Sf) are 
spin operators, ^ is the chemical potential, and the pro¬ 
jection operator P removes zero occupancy in the case of 
the electron doping, i.e., C' la Ci a > 1, while the expo¬ 
nential Peierls factor account for the coupling of electron 
charge to the weak external magnetic field in terms of the 
vector potential A (l). This t-J model © is the strong 
coupling limit of the Hubbard model and the crucial dif¬ 
ficulty of its solution lies in enforcing the local constraint 
of no zero electron occupancy. In the case of the hole 
doping, an intuitively appealing approach to implement 
the local constraint of no double electron occupancy and 
the charge-spin separation scheme is the slave-particle 
formalism^, however, the local constraint of no double 
electron occupancy is explicitly replaced by a global con¬ 
straint in the actual calculations. Following the charge- 
spin separation scheme, a fermion-spin theory has been 
develope d 25 ! 30 , where the local constraint of no double 
electron occupancy can be treated properly in the ac¬ 
tual calculations. To apply the fermion-spin theory to 
the case of the electron doping, the t-J model © can 
be rewritten in terms of a particle-hole transformation 
Cia -A fi_ a as^I, 

h = -tj2 e ~ i[e/h)Ml) 'yU+v°+vY.fU° 

lf)a la 

+ J^TSz-Sz +r p (2) 

If) 

where ff (fia) is the hole creation (annihilation) oper¬ 
ator. This t-J model © in the hole representation is 
subject to the local constraint that double occupancy of 
a site by two fermions of opposite spins is not allowed, 
i.e., 12a fiafi* — 1- The physics of the no double elec¬ 
tron occupancy in the fermion-spin theory is taken into 
account by representing the fermion operator f ia as a 
composite object created b y 25 ' 30 , 

fit = a n s i > /u = a z\^z + j (3) 

where the spinful fermion operator ai a = e _i$,<7 az de¬ 
scribes the charge degree of freedom of the hole together 
with some effects of spin configuration rearrangements 
due to the presence of the doped electron itself (charge 
carrier), while the spin operator Si represents the spin 
degree of freedom of the hole, then the local constraint 
of no double occupancy is satisfied in the actual calcu¬ 
lations. In this fermion-spin representation d3j), the t-J 
model © can be expressed as^i, 

H = t ^e _i ( e / R ) A ( ; h^ ( aJ’ + - t az t S' i + 5zjj_j 5 + a\ + ~^aiiSf S^) 

If) 

M ^ ^ Q j i ( jQ j la H - JeS ^ ^ (4) 

la If) 

where J e ff = (1 — <5) 2 J, and 6 = (aj a ai a ) = (ajai) is the 
charge-carrier doping concentration. At half-filling, the 
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t-J model f4]) is reduced to the AF Heisenberg model on 
a triangular lattice. In the early days of the spin-liquid, it 
was proposed that the strong geometry frustration in the 
triangular-lattice Heisenberg model may completely de¬ 
stroy AFLRQ22. Later, a series of studies with spin-wave 
calculation s 33 ! 34 and numerical simulation !) 35 ’ 36 indicate 
that the triangular-lattice AF Heisenberg model appears 
to have better state with three sublattice magnetic or¬ 
der. However, for the case in the square-lattice, it has 
been shown that AFLRO is destroyed by charge-carrier 
doping with S ~ 0.05 - 0.07 for t/J ~ 2.5 - 522^. It 
is thus possible that the spin-liquid state is attained in 
the triangular-lattice system for sufficiently low doping, 
such as 6 ~ 0.05, due to the strong geometry frustration, 
and then there is no AFLRO in the doped regime where 
superconductivity appears. 

Much of the interest in oxide superconductors is due 
to the fact that these materials represent novel SC mech¬ 
anism for superconductivity. Based on the t-J model 
in the fermion-spin representation, we have developed 
a kinetic-energy-driven SC mechanis m 24 ’ 25 in the case 
without AFLRO for a microscopic description of the SC- 
state of square-lattice cuprate superconductors. This 
kinetic-energy-driven superconductivity is purely elec¬ 
tronic without phonons, where the charge-carrier pairing 
interaction arises directly from the kinetic energy by the 
exchange of spin excitations in the higher powers of the 
doping concentration. In particular, the kinetic-energy- 
driven SC-state is controlled by both the SC gap and 
quasiparticle coherence, then the maximal SC transition 
temperature T c occurs around the optimal doping, and 
decreases in both the underdoped and overdoped regimes. 
On the other hand, since the strong electron correlation 
is common for both these materials^ - — ; these two oxide 
systems may have similar underlying SC mechanism, i.e., 


it is possible that superconductivity in triangular-lattice 
cobaltate superconductors is also driven by the kinetic 
energy. In this case, superconductivity in triangular- 
lattice superconductors has been discussed based on the 
kinetic-energy-driven SC mechanism^!, where although 
the effect from the charge-carrier quasiparticle coherence 
has been dropped, the obtained doping dependence of 
T c is qualitative agreement with experimental daiaXAS, 
of Na^CoC^-y^O. In this paper, we generalize the for¬ 
malism of the kinetic-energy-driven SC mechanism de¬ 
veloped in Ref. dH by considering the effect from the 
charge-carrier quasiparticle coherence, and then apply 
this new form to study the electromagnetic response in 
triangular-lattice superconductors. Following the previ¬ 
ous discussion s 24 ' 25 ’ 31 , the self-consistent equations that 
are satisfied by the full charge-carrier diagonal and off- 
diagonal Green’s functions of the triangular-lattice t-J 
model (HI) in the SC-state at zero magnetic field can be 
obtained explicitly as, 

3(k,w) = 5 {0) (k,w) +g (0) (k,w)[£[ a) (k,w)g(k,w) 

- ^ a) *(k,w)r t (k,w)], (5a) 

r f (k, u>) = 5 {0) (-k, -w)[E^ a) (-k, -w)r t (k,w) 

+ (5b) 

respectively, where gl°) _1 (k, ui) = oj —£k is the mean-field 
(MF) charge-carrier Green’s function, £k = Zt\ 7k — H is 
the MF charge-carrier excitation spectrum, Z is the num¬ 
ber of the nearest-neighbor sites, 7k = X/Z) JA e lk '*> = 

[cos k x + 2 cos(fc x /2) cos(y/3k y /2)]/3, and the spin corre¬ 
lation function \ is defined as \ = , while the 

charge-carrier self-energies have been obtained as, 


J 


E^ a) (k,iw n ) = (Zt) 2 —^ 7 2 +p , +k i^ ff (p + k,ip m + «w n )n(p,p , ,ip m ), (6a) 

p. p / ipm 

E^ a) (k ,iu n ) = (Zt) 2 -^2 ^7 2 +p , +k i^r t (p + k,ip m + «w„)n(p,p , ,ip m ), (6b) 

P P' iPm 


with the spin bubble, 

n(p,p ',ip m ) = 4 ^^£) (0) (p , ,*Pm) 

ip!m 

x £> (0) (p' + pXp’ m + ip m ), (7) 

where D^{1 — l',t — t') = ((S) + (t); Sj7 (f'))) is the MF 
spin Green’s function, and has been evaluated as, 


with the function B p = A [2\ z (ej P — 1) + x{lv ~ e )] > while 
the MF spin excitation spectrum u) p is given by, 

Wp = A 2 ie (Ai - iax z - a+7k^ (e - 7k) 

+ ^2 - J^ ae X ~ aex z 7k^ (1 - £7k) , 

( 9 ) 


D {0) ( P,w) 


B P f 1 

2uj p \ui — ui p 


1 


Cl) + ClV 


with A 1 = aC+(l-a)/(2Z), A 2 = aC z + (1 -a)/(4Z), 
A = 2 ZJ e g, e = 1 + 2 tcj)/J e f[, the charge-carrier’s particle- 
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hole parameter (f> = (a^a^o-), and the spin correlation 
functions X z = {SfS? +fj ), C = (l/Z 2 ) (S+ + -S; + .,), 

and C z = (1 /Z 2 )J2firf{Si+r)Si +r j l )- Since the quantum 
spin operators obey the Pauli algebra, it needs to ap¬ 
ply the decoupling approximation 39,4 - to the higher order 
spin Green’s function for obtaining the MF spin Green’s 
function £)( 0 ^(p,w). In particular, in order to satisfy the 
sum rule of the correlation function (S/~ Sf) = 1/2 in the 
case without AFLRO, an important decoupling parame¬ 
ter a has been introduced in the decoupling approxima¬ 
tion for the higher order spin Green’s function, which can 
be regarded as the vertex correction, and is determined 
self-consistentl y 25 ’ 40 . 
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FIG. 1: (Color online) (a) The spectral intensity maps at the 
charge-carrier Fermi energy in S = 0.20 with T = 0.001J for 
t/J= —2.5. (b) The d-wave gap maps in the Brillouin zone. 


In obtaining Eqs. m and the facts 

E^ a) (-k,-w) = £^, a) (k,w) and rt(-k,-w) = rt(k,w) 

have been used, which indicates that the charge-carrier 
pair gap A^ (uj) = E^ (k, uj) is an even function of uj. 
However, the other charge-carrier self-energy (k, uj) 
is not. It is convenient to break it into its symmetric 
and antisymmetric parts as E^ a \k,w) = E^(k, w) + 
wE[^(k,w), and then £^(k, w) and £^(k, uj) are both 
even function of uj. The antisymmetric part E^ (k, uj) 
of the self-energy E^ a) (k, uj) renormalizes the MF charge- 
carrier spectrum, and is directly related to the charge- 
carrier quasiparticle coherent weight as Z~p(k, uj) = 


1 — Re£^(k, uj), while the symmetric part E^(k, w) of 
the self-energy E^ a ^(k, uj) just renormalizes the chemical 
potential. In this paper, we mainly focus on the low- 
energy behavior, and in this case, the charge-carrier pair 
gap and quasiparticle coherent weight can be generally 
discussed in the static limit, i.e., = E^^k, w) |^=o, 

and Z~p (k) = 1 - ReE^ a) (k,w) | w=0 . Although Z aF (k) 
still is a function of momentum, however, as a qualitative 
discussion, the wave vector k in Z aF ( k) can be chosen in 
the high symmetry point of the Brillouin zone, i.e., 

-J— = 1 ~ ReE^(k,cu = 0) | ko , (10) 

^aF 

with ko = [7 t/ 3, \/37r/3]. In triangular-lattice super¬ 
conductors, a central issue is whether the charge-carrier 
pair gap has nodes at the Fermi surface or not. Ex¬ 
perimentally, it is far from reaching a consensus on the 
pairing symmetry. From early NMR and NQR mea¬ 
surements, the contradictory results were obtained, since 
some experimental data are consistent with the case of 
the existence of a coherence peak indicating a complete 
gap over the Fermi surfac o 10 ! 11 , while other experimen¬ 
tal results suggest no coherence pea k — i— . In particu¬ 
lar, although the recent experimental results— obtained 
from the specific-heat measurements do not give unam¬ 
biguous evidence for either the presence or absence of 
the nodes in the energy gap, the experimental data of 
the specific-heal— are consistent with these fitted re¬ 
sults obtained from phenomenological Bardeen-Cooper- 
Schrieffer (BCS) formalism with the d-wave (di+fc ?2 pair¬ 
ing) symmetry without gap nodes. On the theoretical 
hand, according to the irreducible representations of the 
triangular-lattice system, it has been pointed out that 
there are three possible basis functions of even parity^, 
i.e., one s-like function Sk = cosk x + cos[(k x — \/3k y )/2] + 
cos[(fc x + y/3k y )/2], and two d-like functions, dik = 
2cosk x —cos[(k x —\/3ky)/2]—cos[(k x +V3ky)/2] andcfck = 
V / 3cos[(fc a; + v / 3fc y )/2] — \/3cos[(fc x — y/3k y )/2). However, 
with the different linear combinations of these basis func¬ 
tions, one find— based on the variational Monte Carlo 
simulation of the resonating-valence-bond wave function 
that the lowest energy state of the AF triangular-lattice 
Heisenberg model is the d-wave (dik + ic? 2 k) state with 
the energy gap Ak oc A(dik + icfok)- Furthermore, it 
has been shown based on the numerical simulations that 
this d-wave state also is the lowest state around the 
electron-doped regime where superconductivity appears 
in triangular-lattice superconductors^—. In particular, 
the recent theoretical studies based on a phenomenolog¬ 
ical analysis^ and a combined cluster calculation and 
renormalization group approach— show that this d-wave 
state naturally explains some SC-state properties as in¬ 
dicated by experiments. In this case, we only consider 
the case with the d-wave pairing symmetry, 

A^=A W(d lk + id 2 k), (11) 

and then the full charge-carrier diagonal and off-diagonal 
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Green’s functions in Eq. © can be obtained explicitly 
as, 


g{ k,w) 
r f (k, w) 


Z a F 


~Z aF 


II 2 

U aU. 


Kk 


^ -G i k & T A lk , 

^( a ) / i 

—‘^l. / -L 


A Zk 


(12a) 
,(12b) 


with the charge-carrier quasiparticle energy spectrum 
E a k = y£ k + | | 2 , the renormalized charge-carrier 

excitation spectrum £ k = Z a p£k, and the renormalized 
charge-carrier pair gap = -ZapA^ 1 , while the charge- 
carrier quasiparticle coherence factors, 


uL = 


H 1+ li 


(13a) 

(13b) 


satisfy the constraint U 2 k + V 2 k = 1 for any wave vec¬ 
tor k. In spite of the pairing mechanism driven by the 
kinetic energy by the exchange of spin excitations, the 


J 


result in Eq. in is a standard BCS expression for a 
charge-carrier d-wave pair state. It should be empha¬ 
sized that in triangular-lattice superconductors, there is 
a large charge-carrier Fermi surface around the T point in 
the Brillouin zone as well as six small hole-pockets near 
the K points as shown in Fig. [T}i. In particular, the nodes 
of the charge-carrier d-wave pair gap HID exist only on 
the six hole-pockets and not on the large charge-carrier 
Fermi surface as shown in Fig. Qj>. However, these nodes 
around the six hole-pockets are far from the large charge- 
carrier Fermi surface, and therefore there are no gapless 
charge-carrier quasiparticle excitations. In this case, the 
effect from these charge-carrier quasiparticles around the 
nodes is unimportant, since everything happens at the 
charge-carrier Fermi surface. It is thus expected that 
the basic behavior of the evolution of the magnetic-field- 
penetration depth (then the superfluid density) with tem¬ 
perature in triangular-lattice superconductors is much 
different from that in square-lattice superconductors. 

According to the full charge-carrier Green’s functions 
in Eq. (1121) and spin Green’s function in Eq. ©, now 

the charge-carrier self-energies E^(k, w) and E^k, u>) 
can be evaluated explicitly as, 


E^(k, w ) = 


]\f2 “pp'k 


pp'n 


II 2 

^ap+k 


UJ 


V 2 

*ap+k 


E 


,(n) 


lapp'k 


UJ — UJ. 


J npp' 


- F.„ 


^ap+k 

pp'n ^ 

A n ) 


p(n) 

lapp'k 


F (n) 

± 2app'k 


^npp' -£^ap+k ^ ^npp' -^ap+k 


R 


lapp'k 


UJ — UJr 



R 


,(n) 


2app'k 


UJ ~ UJ. 


npp 


— E, 


'ap+k 


(14a) 


(14b) 


respectively, with n = 1,2, flpp, k = Z a F(Ztj p+p > +k ) 2 B P ' B P+P J (Auj p >lo p+p >), u npp ' = w p+p - - (-l) n Wp-, and the 
functions, 


F ia ] pp'k = «F(^ap+k){l + n B (w p ' +p ) + n B [(-l)" +1 w p /]} + n B (cv+p W(-l) n+1 ^p'], (15a) 

= [ 1-n F( £: ap+k)]{l + nB(wp/ + p)-|-n B [(—l)" +1 cjp/]} + n B (wp/ +p )n B [(-l)" +1 cjp/], (15b) 

where n B (w) and n F (uj) are the boson and fermion distribution functions, respectively. In this case, the charge-carrier 
quasiparticle coherent weight Z aF in Eq. (flOl) and charge-carrier pair gap parameter in Eq. (fill) satisfy following 
two self-consistent equations, 


1 

Z aF 

1 


N 2 


E(-i)" +, A‘A. 


F 


( n) 

lapp'ko 


T?( n ) 

2app'k 0 


pp'n 


(a^npp' Akp+ko ) 2 (cUnpp' T A a p-|-ko 'j 2 


r A ( d ) * A ( d ) 

_ ('-If n 7 o( a ) k A p+k 

N3 R ( ^ ^ aF pp'k Ea 

pp'kn 


E7’( n ) 

lapp'k 


p(«) 

‘ 2app'k 


^ap+k 


. ^npp' ^ap+k 


^npp' H - ^ap+k 


(16a) 

(16b) 


respectively, with A 


(d) _ 

k — 


dik + *d 2 k- These two equations 


(I16al) and (I16bl) must be solved simultaneously with fol- 
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lowing self-consistent equations, 

0 = i 1 ~ ,(17a) 

k x 7 

s = wS z -( 1 -^: , “ h i5' 5E -i)' (17b) 

k x 7 


* = < 17c > 

k 

c = ( 17d > 

k 

l = (i7e) 

k 

^ = ^E^|f^coth[i/3 Wzk ], (m) 

k 

c ' = (17g) 

k 


then all the order parameters, the decoupling parame¬ 
ter a, and the chemical potential /i are determined self- 
consistently without using any adjustable parameters. 



FIG. 2: The charge-carrier pair gap parameter as a function 
of doping with T = 0.001 J for t/J = —2.5. 

These equations in Eqs. © and © have been calcu¬ 
lated self-consistently, and the result of the charge-carrier 
pair gap parameter A( a ) as a function of doping for pa¬ 
rameter t/J = —2.5 with temperature T = 0.001J is 
shown in Fig. [2] It is shown clearly that the charge- 
carrier pair gap parameter takes a domelike shape 
with the underdoped and overdoped regimes on each side 
of the optimal doping <5 op timai ~ 0.19, where A^ a ) reaches 
its maximum. Moreover, we have made a series of calcu¬ 
lations for A< a ) at different temperatures, and the result 
shows that the charge-carrier pair gap parameter A( a ) fol¬ 
lows qualitatively a BCS-type temperature dependence. 

T c on the other hand can be obtained self-consistently 
from the self-consistent equations © and © by 



FIG. 3: T c as a function of doping for t/J = —2.5. Inset: the 
corresponding experimental result of Naa;Co 02 -yH 20 taken 
from Ref. @ 

the condition A (a ) = 0, and the result of T c as 
a function of doping for t/J = —2.5 is plotted in 
Fig. [3] in comparison with the corresponding experi¬ 
mental result^ of Na^CoC^-y^O (inset). Obviously, 
the experimental result^ 10 of the doping dependence 
of T c in the triangular-lattice cobaltate superconductor 
Na^CoCVy^O is qualitatively reproduced. The opti¬ 
mal T c occurs in a narrow range of doping, and then de¬ 
creases for both underdoped and overdoped regimes, in 
dramatic analogy to the phase diagram of square-lattice 
cuprate superconductors^. However, in comparison with 
the result of T c obtained from square-lattice cuprate su¬ 
perconductors, the present result also shows that the ge¬ 
ometric frustration, accompanied by large fluctuations, 
suppresses T c to low temperatures. 

Within the above framework of the kinetic-energy- 
driven superconductivity, we now turn to discuss the dop¬ 
ing and temperature dependence of the electromagnetic 
responses in triangular-lattice superconductors. For dis¬ 
cussions of the Meissner effect in a superconductor, one 
usually starts from the general relation between the cur¬ 
rent and the vector potentia l 17 ’ 49 , 

J M (q,w) = -'Y^K, lv (q,w)A v {q,u) (18) 

V 

where the Greek indices label the axes of the Carte¬ 
sian coordinate system, while the nonlocal kernel of 
the response function K^(q,oj) can be expressed as 
K^ v { q,w) = K$(q,u) + K$(q,u), with K$(q,uj) 
and Kj%)(q, u>) that are the corresponding diamagnetic 
and paramagnetic parts. 

In the fermion-spin representation d3j) , the vector po¬ 
tential A has been coupled to the electron charge which 
are now represented by /q = and /q = aqS) + - In 

this case, the electron polarization operator is expressed 
as P = e ^l^ierCia = e '/2i Rt ®; a U and then the cor¬ 
responding electron current operator is obtained by eval¬ 
uating the time derivative of this polarization operator 
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j = dP/dt = i[H, P]/ft2S. In the linear response ap¬ 
proximation with respect to A U {1), this electron current 
operator is reduced as j = j (d ) +j^ p - ) , with the correspond¬ 
ing diamagnetic (d) and paramagnetic (p) components of 
the electron current operator are given by, 

j (d) = 

It) 

+ aj+fnaiiS j-, (19a) 

j (p) = 

Irj 

+ a ]+f)i a n s r s tv) > ( 19b ) 

respectively. The diamagnetic component of the electron 
current operator in Eq. (1 19 all is proportional to the vec¬ 
tor potential, and therefore the diamagnetic part of the 
response kernel can be obtained directly as, 

6e 2 1 

K$( q,w) = --^ 2 (20) 

with the London penetration depth A^ 2 = — 6e 2 \</>i//i 2 . 

The paramagnetic part of the response kernel, on 
the other hand, is directly related to the elec¬ 
tron current-current correlation function P yv (q, r) = 

— (T T {jj?\q, T)jl P \— q, 0)}), and can be expressed as 
= P y „(q,uj). In the fermion-spin approach, 
the paramagnetic component of the electron current op¬ 
erator in Eq. (Il9bl) can be decoupled as, 

.(p) i e Xt - t 

y p> = — j-2^W J l+fl<J a la 

lr)(j 

ifj 


As in the case of square-lattice superconductors^, the 
second term in the right-hand side refer to the contri¬ 
bution from the electron spin, and can be shown that 
Sf + ~+Sf -) = 0, i.e., there is no direct contri- 

lr) 

bution for the electron current-current correlation func¬ 
tion P At ^(q, r) from the electron spin, and then the ma¬ 
jority contribution for P /Jl/ (q, r) comes from the electron 
charge, however the strong interplay between charge car¬ 
riers and spins has been considered through the spin’s or¬ 
der parameters entering in the charge-carrier part of the 
contribution to the current-current correlation P M „(q, r). 

The density operator is summed over the position 
of all particles, i.e, p{l) = -e[l/N] C} a Cio = 
—e[l/(2 N)\Y^i a a\,j a hy, and then its Fourier transform 
can be expressed as p(q) = —e/(2 N) ]P k(T (i^flk+qa- For 
a convenience in the following discussions, the paramag¬ 
netic component of the electron current operator in Eq. 
(Il9bl) and density operator can be rewritten into the four- 
current operator in the Narnbu representation in terms 
of the charge-carrier Nambu operators Tj) = (a\^,a-k\) 

and ^ k+q = (a fc+gt ,a 1 ) fc _ l?i ) T as, 


4 P) (q) = kT7i(k,k + q)tf k+q , (22) 

k 


(21) with the bare current vertex, 

I 


7 /i (k, k -I- q) = < 


-l^ei* qx {cosjq x sm(k x +lq x )+cos&q y sm(±k x +jq x )cos(&ky+&q y ) 

+i[sm±q x sm(k x +±q x )+sm^q y c 0 s(±k x +lq x )sm(^k y +^q y )]} 

[ cos iq xCOS (ik x +jq x )sm(^ky+^q y ) 

+ismjq x sm(^k x +jq x )cos(Y-k y +^fq y )] 

-| er 3 


for p=x 

(23) 

for p=y 
for //=0 


In this case, the current-current correlation function is obtained as, 


P»v( q, w„) = ^2 7/i( k + q, k)y*(k + q, k)i ^ Tr[<j(k + q ,iw n + iu m )g(li,iu m )}, 


(24) 


where the full charge-carrier Green’s function g(k,u>) in the Nambu representation can be expressed in terms of the 
full charge-carrier Green’s function m as, 


<?( k,w) 


ZaF 


WT 0 + £ k T 3 


- A^ a) (di k ri + d 2 ^T 2 ) 


(25) 













Substituting this charge-carrier Green’s function (l25l) into Eq. (12411 . the paramagnetic part of the response kernel in 
the static limit (w ~ 0) is evaluated as, 


^(q,0) = — ^ 7 M (k + q,k) 7 *(k + q,k)[Li a) (k,q)+4 a) (k,q)] = K$(q,0)5^, 

k 

with the functions l['^ (k, q) and (k, q) are given by, 

L^(k,q) = Z 2 


(26) 


/ c c iaWaW* i i*W*aW ' 

( a )n, _ 7 2 If sk?k+q + 2 ZA Zk ZA Zk+q + 2 ^Zk ^Zk+q 


nF(E/ a k) ^F(E/ a k+q) 


J aF 


4 a) (k,q) = Z, 


^'ak-S'ak+q 
1 A ( a ) A ( a )* _L 1 A ( a )* A ( a ) 


-^ak -S'ak+q 


2 ^k+q+ 2A Z kAzk+q+ 2 A Zk A Zk+q\ «F(-^ak) + «F(A a k+q) — 1 

aF I 77> I TT \ TT 7 


^'ak-^'ak+q 


-^ak -^ak+q 


(27a) 

(27b) 


respectively. In this case, the kernel of the response function in Eq. m is now obtained from Eqs. and m as, 


A V(q>°) = 


^+^P)(q,0) 


l_' v L 


<W 


(28) 


In the long-wavelength limit, i.e., |q| —> 0, the function [k, q —> 0) vanishes, then the paramagnetic part of the 
response kernel in Eq. m is reduced as, 


tf<P>(q->0,0) = Z; 


24e 2 1 

~Wn 


Y.xV™ 2 ( l Msm\fk y ) lL mw 

, z Z q - ^0 -C^ak -^ak+q 


(29) 


However, at zero temperature (T = 0), A'yy^(q —> 0,0)|t=o = 0, and then the long-wavelength electromagnetic 
response is determined by diamagnetic part of the response kernel only. On the other hand, at T = T c , the 
charge-carrier gap parameter A^\t=t c = 0, and in this case, the paramagnetic part of the response kernel in the 
long-wavelength limit can be evaluated as, 


K &(q ^ °’ 1 °) = Z aF^^ E X 2 * 2 cos2 (^*) si * 2 (^k y ) lim nF( ^_^ k+q) 


A V 


(30) 


which exactly cancels the diamagnetic part of the response kernel in Eq. (l20l) . and then the Meissner effect in 
triangular-lattice superconductors disappears for all temperatures T >T C . These results also reflect that the Meissner 
effect is strongly temperature dependent. To show this point clearly, we introduce an effective superfluid density 
n s (T) at temperature T, which is defined in terms of the paramagnetic part of the response kernel as, 


4 p , } (q^o,o) 




n s {T) 
n s ( 0) 


(31) 


where the ratio n s (T)/n s (0) of the effective superfluid densities at temperature T and zero-temperature is obtained 
directly from the paramagnetic part of the response kernel (1291) as, 


n s {T) 
n s ( 0) 


1 - a(Af 


24e 2 1 
~h?~N 



(g/3£ak 1)2 


(32) 


In this case, the kernel of the response function in Eq. 
(E51) can be expressed explicitly in terms of the effective 
superfluid density as, 

AV(q^0,0) = (33) 

a l n a {\J) 

In Fig. 0] we plot this effective superfluid density 


n s (T)/n s ( 0) as a function of temperature at 6 = 0.15 
for t/J = —2.5, where n s (T)/n s ( 0) decreases with in¬ 
creasing temperatures, and vanishes at T c , then all the 
charge-carriers are in the normal fluid for temperatures 
T > T c . To sum up, within the kinetic-energy driven 
SC mechanism, we find that: (a) the Meissner effect 
in triangular-lattice superconductors is obtained for all 
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temperatures T < T c throughout the SC dome; (b) the 
electromagnetic response kernel goes to the London form 
in the long-wavelength limit [see, e.g., Eq. ([33)l ]: (c) al¬ 
though the electromagnetic response kernel (1281) is not 
manifestly gauge invariant within the bare current vertex 
(123p . however, we can keep the gauge invariance within 
the dressed current vertex as it has been done in the case 
for square-lattice superconductors^; (d) in spite of the 
pairing mechanism driven by the kinetic energy by the ex¬ 
change of spin excitations, the kinetic-energy-driven SC- 
state in triangular-lattice superconductors still is conven¬ 
tional BCS-like with the d-wave symmetry without gap 
nodes at the charge-carrier Fermi surface, which leads to 
a fact that the superfluid density therefore follows essen¬ 
tially a BCS-type temperature dependence as in the case 
of conventional superconductors^. 



FIG. 4: The effective superfluid density as a function of tem¬ 
perature at 5 = 0.15 for t/J = —2.5. 


III. DOPING DEPENDENCE OF 
ELECTROMAGNETIC RESPONSE 

In the above discussions, it is shown that once the re¬ 
sponse kernel K^ is known, the effect of a weak elec¬ 
tromagnetic field on a superconductor can be quantita¬ 
tively characterized by experimentally measurable quan¬ 
tities such as the magnetic-field-penetration depth and 
superfluid density. However, the result of the effective 
superfluid density in Eq. (15^1) obtained from the response 
kernel in Eq. (1331) can not be used for a direct comparison 
with the corresponding experimental data of triangular- 
lattice superconductors because the kernel function de¬ 
rived within the linear response theory describes the re¬ 
sponse of an infinite syste m 23 ' 50 , whereas in the actual 
problem of the penetration of the field and the system 
has a surface, i.e., it occupies a half-space x > 0. In 
this case, we need to impose boundary conditions for 
charge carriers, which can be done within the simplest 
specular reflection mode&S 2 with a two-dimensional ge¬ 
ometry of the SC plane. Following our previous discus¬ 
sions of the electromagnetic response in square-lattice 


superconductor s 23 ' 50 , the local magnetic field profile of 
triangular-lattice superconductors can be evaluated ex¬ 
plicitly as, 


h z (x) = — 

7r 


dq a 


q x sin (q x x) 


^oK yy (q x , 0, 0) + q x 


(34) 


and then the magnetic-field-penetration depth is ob¬ 
tained from this local-magnetic-field profile as, 

1 r°° 

A (T) = — / h z (x)dx 

2 
7r 

For a convenience in the following discussions, we intro¬ 
duce a characteristic length scale a o = yWa/jjLo^J . Us¬ 
ing the lattice parameter a ss 0.282 nrn for Na^CoC^- 
yH 2 0, this characteristic length is obtained as ao ~ 83.9 
nm. 


d q x 


1-loKyy (q x , 0, 0) + q x 


(35) 



FIG. 5: The magnetic-field-penetration as a function of tem¬ 
perature at 5 = 0.15 for t/J = —2.5. Inset: the magnetic- 
field-penetration difference AA(T) = A(T) — A(0) (solid line) 
as a function of temperature at <5 = 0.15 for t/J = —2.5, while 
the dashed line is obtained from a numerical fit AA (T) = 
Aexp[-BA (a \T)/T] with A ~ 2789.62 and B ~ 0.59. 


We are now ready to discuss the electromagnetic re¬ 
sponse in triangular-lattice superconductors. At zero 
temperature, the obtained magnetic-field-penetration 
depths from Eq. (l35l) are A(0) ~ 362.96 nm, A(0) ~ 
316.38 nm, and A(0) « 294.36 nm for 6 = 0.15, 6 = 0.17, 
and <5 = 0.19, respectively. However, at T = T c , since the 
kernel of the response function K ^ v {q —► 0, 0)| t=t c = 0, 
the magnetic-field-penetration depth is found as A(T c ) = 
oo, i.e., the external magnetic field can penetration 
through all the main body of the system for T >T C , and 
then the Meissner effect does not exist in the normal- 
state. On the other hand, A(T) is sensitive to low-lying 
excitations. To show this point clearly, A(T)/A(0) as a 
function of temperature at S = 0.15 for t/J = —2.5 is 
plotted in Fig. [5] It is seen that below temperatures 
T < 0.25T c , A(T) is practically independent tempera¬ 
ture, which is a reflection of the absence of the d-wave gap 
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nodes at the large charge-carrier Fermi surface. However, 
above temperatures T > 0.25T C , A(T) increases rapidly 
with increasing temperature. In particular, we have fit¬ 
ted our present theoretical result of the magnetic-field- 
penetration depth difference AA (T) = A(T) — A(0), and 
the fitted result is shown in inset of Fig. [5] We thus find 
that AA (T) vary exponentially as a function of temper¬ 
ature (AA (T) = Aexp[-BA^(T)/T] with A ~ 2789.62 
and B ~ 0.59), which is expected result in the case 
without the d-wave gap nodes at the large charge-carrier 
Fermi surface. 

An external magnetic field acts on the SC-state of 
triangular-lattice superconductors as a perturbation. In 
the linear response form (1181) . the nonlocal relation be¬ 
tween the supercurrent and the vector potential in the 
coordinate space holds due to the finite size of charge- 
carrier pairs. In particular, the size of charge-carrier 
pairs is of the order of the coherence length £(k) = 
where up = fi _1 <9£k/dk|fc p is the charge- 
carrier velocity at the large charge-carrier Fermi surface, 
which shows that the size of charge-carrier pairs is mo¬ 
mentum dependent. In general case, although the exter¬ 
nal magnetic field decays on the scale of the magnetic- 
field-penetration length A(T), any nonlocal contributions 
to measurable quantities are of the order of k~ 2 , where 
the Ginzburg-Landau parameter n is the ratio of the 
magnetic-field-penetration depth A and the coherence 
length (. However, for the charge-carrier d-wave pair 
gap m, there is no gap nodes at the large charge-carrier 
Fermi surface. In this case, the momentum dependent co¬ 
herence length £(k) can be replaced approximately with 
the isotropic one £ 0 = Hvf/(ttA^), and then the con¬ 
dition for the local limit is satisfied. As a consequence, 
triangular-lattice superconductors are type-II supercon¬ 
ductors due to the existence of the anisotropic energy 
gap over the large charge-carrier Fermi surface, where 
nonlocal effects are negligible, and then the electrody¬ 
namics is purely local and the magnetic field decays ex¬ 
ponentially over a length of the order of a few hundreds 
nm. In this local limit, the pure d-wave pairing state in 
the kinetic-energy-driven SC mechanism gives a tempera¬ 
ture dependence of the magnetic-field-penetration depth 
as AA(T) oc exp[—A^ a ^(T)/T]. This is much different 
from the case in square-lattice cuprate superconductors, 
where the characteristic feature of the d-wave charge- 
carrier pair gap is the existence of the four nodes at 
the charge-carrier Fermi surface, and then the quasi¬ 
particle excitations are gapless and affect particularly 
the physical properties at the extremely low tempera¬ 
tures. These gapless quasiparticle excitations in square- 
lattice cuprate superconductors lead to a divergence of 
the coherence length £(k) around the gap nodes, and 
then the behavior of the temperature dependence of the 
magnetic-field-penetration depth depends sensitively on 
the quasiparticle scattering. At the extremely low tem¬ 
peratures, the quasiparticles selectively locate around the 
gap nodal region, and then the major contribution to 
measurable quantities comes from these quasiparticles. 


In this case, the Ginzburg- Landau ratio ft(k) around the 
gap nodal region is no longer large enough for the sys¬ 
tem to belong to the class of type-II superconductors, 
and the condition of the local limit is not satisfie d ! !— , 
which leads to the system in the extreme nonlocal limit, 
and therefore the nonlinear behavior in the tempera¬ 
ture dependence of the magnetic-field-penetration depth 
is observed experimentally!!. On the other hand, with 
increasing temperatures, the quasiparticles around the 
gap nodal region become excited out of the condensate, 
and then the nonlocal effect fades away, which leads 
to that the magnetic-field-penetration depth crossovers 
to the linear temperature dependence!!. However, the 
present result of the temperature dependence of the 
magnetic-field-penetration depth in triangular-lattice su¬ 
perconductors is very similar to the case in conventional 
superconductors!!, where the characteristic feature is the 
existence of the isotropic energy gap at the Fermi surface, 
and then the temperature dependence of the magnetic- 
field-penetration depth exhibits an exponential behavior. 



FIG. 6: The superfluid density as a function of doping with 
T = 0.001J for t/J = -2.5. 

Now we turn to discuss the doping dependence of the 
superfluid density p s (T ), which is a measure of the phase 
stiffness, and is defined in terms of the magnetic-field- 
penetration depth A(T) as p s (T) = 1/A 2 (T). In this case, 
we have performed firstly a calculation for the doping de¬ 
pendence of p s in triangular-lattice superconductors for 
all levels of doping throughout the SC dome, and the 
result is plotted in Fig. [6] In analogy to the dome¬ 
like shape of the doping dependence of T c shown in Fig. 
[3j p s also displays a dome-like shape of the doping de¬ 
pendence, i.e., it increases with increasing doping in the 
lower doped regime, and reaches a maximum (a peak) 
around the critical doping Critical ~ 0.21, then decreases 
in the higher doped regime. Moreover, p s of triangular- 
lattice superconductors in the underdoped regime van¬ 
ishes more or less linearly with decrease of the charge- 
carrier doping concentration 6. In square-lattice cuprate 
superconductors, one of the most unconventional natures 
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is that T c scales with p s following the so-called Uemura 
relation as T c oc const x p s in the underdoped regime^. 
It is interesting to know if triangular-lattice supercon¬ 
ductors also obey this relation. In this case, we have 
fitted the relation between T c and p s in the underdoped 
regime, and the result shows that triangular-lattice su¬ 
perconductors satisfy the similar Uemura relation in the 
underdoped regirn o 26 ' 27 . Incorporating the result ob¬ 
tained from square-lattice cuprate superconductors^, it 
thus implies that the Uemura relation may be a universal 
relation in strongly correlated superconductors in despite 
of whether the gap nodes at the charge-carrier Fermi sur¬ 
face exist or not. 

The essential physics of the dome-like shape of the 
doping dependence of p s in triangular-lattice supercon¬ 
ductors is the same as in the case of square-lattice 
superconductor s 23 ' 25 , and also can be attributed to the 
dome-like shape of the doping dependence of A^ a ^. This 
follows a fact that p s (T ) in triangular-lattice supercon¬ 
ductors is intriguingly related to the current-current cor¬ 
relation function, and therefore the variation of the su¬ 
perfluid density with doping is coupled to the doping de¬ 
pendence of the charge-carrier pair gap parameter A^ a b 
In particular, the charge-carrier pair gap parameter A^ a ^ 
measures the strength of the binding of two charge car¬ 
riers into a charge-carrier pair. On the other hand, the 
superfluid density p s is a measurement of the phase stiff¬ 
ness, and is proportional to the squared amplitude of 
the charge-carrier pair macroscopic wave functions. In 
this case, both p s and A^ describe the different aspects 
of the same charge-carrier quasiparticles, and then the 
dome-like shape of the doping dependence of p s in Fig. 
El is a natural consequence of the dome-like shape of the 
doping dependence of A( a ) shown in Fig. [2] 



FIG. 7: The superfluid density as a function of tempera¬ 
ture at S = 0.15 for t/J = —2.5. The experimental result 
of Na^CoCU • yFhO (solid squares) taken from Ref. [53. 

The superfluid density shown in Fig. [G] also is 
strongly temperature dependence. When the temper¬ 
ature T = T c , the kernel of the response function 
—> 0,0)|t=t c = 0, and then A(T c ) = oo as men¬ 


tioned above, which leads to p s {T c ) = 0, and is consis¬ 
tent with the result of the effective superfluid density 
obtained from Eq. (l32ll . For a better understanding 
of the basic behavior of p s (T) as a function of temper¬ 
ature, we have made a series of calculations for p s (T) 
at different temperatures, and the result of p s {T) as a 
function of temperature at S = 0.15 for —t/J = 2.5 is 
plotted in Fig. [7] in comparison with the corresponding 
experimental result^! of Na x Co0 2 • j/H 2 0 (solid square). 
Our present calculations thus qualitatively reproduce the 
overall evolution of the superfluid density with tempera¬ 
ture in Na^CoCU • yH 2 Q2£. In corresponding to the re¬ 
sult of the temperature dependence of the magnetic-field- 
penetration depth shown in Fig. [5l p s (T) is also indepen¬ 
dence of temperature below temperatures T < 0.25T c , 
and then decreases dramatically with increasing temper¬ 
ature for temperatures T > 0.25T C , eventually vanishing 
together with superconductivity at T c . The calculation 
based on the kinetic-energy-driven SC mechanism with 
the d-wave charge-carrier pair gap CD thus gives a good 
agreement with the observed superfluid density data of 
Na x Co0 2 • j/H 2 0. 

IV. CONCLUSIONS 

Within the framework of the kinetic-energy driven su¬ 
perconductivity, we have performed a calculation of the 
doping and temperature dependence of the Meissner ef¬ 
fect in triangular-lattice superconductors for all temper¬ 
atures T < T c throughout the SC dome. Our results 
indicate that the magnetic-field-penetration depth shows 
an exponential temperature dependence due to the ab¬ 
sence of the d-wave gap nodes at the large charge carrier 
Fermi surface. In particular, the experimental result of 
the temperature dependence of the superfluid density in 
cobaltate superconductors can be qualitatively described 
in terms of the d-wave pairing state. However, as a nat¬ 
ural consequence of the dome-like shape of the doping 
dependence of the charge-carrier pair gap parameter and 
T c , the superfluid density increases with increasing dop¬ 
ing in the lower doped regime, and reaches a highest value 
(a peak) around the critical doping, then decreases in the 
higher doped regime. 

Acknowledgments 

This work was supported by the funds from the Min¬ 
istry of Science and Technology of China under Grant 
Nos. 2011CB921700 and 2012CB821403, and the Na¬ 
tional Natural Science Foundation of China under Grant 
Nos. 11274044 and 11447144. JQ is supported by the 
Fundamental Research Funds for the Central Universi¬ 
ties under Grant No. FRF-TP-14-074A2, and the Bei¬ 
jing Higher Education Young Elite Teacher Project under 
Grant No. 0389. 







12 


E-mail: spfeng@bnu.edu.cn 

2 See, e.g., the review, Tom Timusk and Bryan Statt, Rep. 
Prog. Phys. 62 , 61 (1999). 

3 See, e.g., the review, Hideo Aoki, J. Phys.: Condens. Mat¬ 
ter 16 , VI (2004). 

4 J. G. Bednorz and K. A. Muller, Z. Phys. B 64, 189 (1986). 

5 A. P. Ramirez, R. J. Cava, J. J. Krajewski, and W. F. Peck, 
Jr., Phys. Rev. B 49 , 16082 (1994); R.J. Cava, H.W. Zand- 
bergen, A.P. Ramirez, H. Takagi, C.T. Chen, J.J. Krajew¬ 
ski, W.F. Peck Jr., J.V. Waszczak, G. Meigs, R.S. Roth, 
L.F. Schneemeyer, J. Solid State Chem. 104 , 437 (1993). 

6 Kazunori Takada, Hiroya Sakurai, Eiji Takayama- 
Muromachi, Fujio Izumi, Ruben A. Dilanian, and 
Takayoshi Sasaki, Nature 422 , 53 (2003). 

1 R. E. Schaak, T. Klimczuk, M. L. Foo, and R. J. Cava, 
Nature 424 , 527 (2003). 

8 C. J. Milne, D. N. Argyriou, A. Chemseddine, N. Aliouane, 
J. Veira, S. Landsgesell, and D. Alber, Phys. Rev. Lett. 93, 
247007 (2004). 

9 Hiroya Sakurai, Naohito Tsujii, Osamu Suzuki, Hideaki 
Kitazawa, Giyuu Kido, Kazunori Takada, Takayoshi 
Sasaki, and Eiji Takayama-Muromachi, Phys. Rev. B 74 , 
092502 (2006). 

10 C. Michioka, H. Ohta, Y. Itoh, and K. Yoshimura, J. Phys. 
Soc. Jpn. 75 , 063701 (2006). 

11 Yoshiaki Kobayashi, Mai Yokoi, Masatoshi Sato, J. Phys. 
Soc Jpn. 72 , 2453 (2003). 

12 Tatsuya Fujimoto, Guo-qing Zheng, Y. Kitaoka, R. L. 
Meng, J. Cmaidalka, and C. W. Chu, Phys. Rev. Lett. 
92 , 047004 (2004). 

13 Yoshihiko Ihara, Kenji Isliida, Hideo Takeya, Chishiro Mi¬ 
chioka, Masaki Kato, Yutaka Itoh, Kazuyoshi Yoshimura, 
Kazunori Takada, Takayoshi Sasaki, Hiroya Sakurai, and 
Eiji Takayama-Muromachi, J. Phys. Soc. Jpn. 75 , 013708 
(2006). 

14 M. Kato, C. Michioka, T. Waki, Y. Itoh, K. Yoshimura, K. 
Ishida, H. Sakurai, E. Takayama-Muromachi, K. Takada, 
and T. Sasaki, J. Phys.: Condens. Matter 18, 669 (2006). 

15 Guo-qing Zheng, Kazuaki Matano, R. L. Meng, J. Cmaid¬ 
alka, and C. W. Chu, J. Phys.: Condens. Matter 18 , L63 
(2006). 

16 Oleg A. Starykh, Rep. Prog. Phys. 78 , 052502 (2015). 

17 J. R. Schrieffer, Theory of Superconductivity, Benjamin, 
New York, 1964. 

18 See, e.g., the review, B. A. Bonn and W. N. Hardy, in 
Physical Properties of High Temperature Superconductors 
V, edited by D. M. Ginsberg (World Scientific, Singapore, 
1996), p. 7. 

19 See, e.g., the review, C. C. Tsuei and J. R. Kirtley, Rev. 
Mod. Phys. 72 , 969 (2000). 

20 Y. J. Uemura, L. P. Le, G. M. Luke, B. J. Sternlieb, W. 
D. Wu, J. H. Brewer, T. M. Riseman, C. L. Seaman, M. 
B. Maple, M. Ishikawa, D. G. Hinks, J. D. Jorgensen, G. 
Saito, and H. Yamochi, Phys. Rev. Lett. 66, 2665 (1991). 

21 A. Suter, E. Morenzoni R. Khasanov, H. Luetkens, T. 
Prokscha, and N. Garifianov, Phys. Rev. Lett. 92 , 087001 
(2004); S. Kamal, Ruixing Liang, A. Hosseini, D. A. Bonn, 
and W. N. Hardy, Phys. Rev. B. 58 , R8933 (1998); Shih- 
Fu Lee, D. C. Morgan, R. J. Ormeno, D. M. Broun, R. A. 
Doyle, J. R. Waldram, and K. Kadowaki, Phys. Rev. Lett. 
77 , 735 (1996); D. M. Broun, W. A. Huttema, P. J. Turner, 


S. Ozcan, B. Morgan, Ruixing Liang, W. N. Hardy, and D. 
A. Bonn, Phys. Rev. Lett. 99 , 237003 (2007); Mun-Seog 
Kim, John A. Skinta, Thomas R. Lemberger, A. Tsukada, 
and M. Naito, Phys. Rev. Lett. 91 , 087001 (2003); T. R. 
Lemberger, I. Hetel, A. Tsukada, M. Naito, and M. Ran- 
deria, Phys. Rev. B 83 , 140507(R) (2011). 

22 S. K. Yip and J. A. Sauls, Phys. Rev. Lett. 69 , 2264 (1992); 
M. Franz, I. Affleck, and M. H. S. Amin, Phys. Rev. Lett. 
79 , 1555 (1997); M. H. S. Amin, M. Franz, and I. Af¬ 
fleck, Phys. Rev. Lett. 84 , 5864 (2000); Mei-Rong Li, P. 
J. Hirschfeld, and P. Woffle, Phys. Rev. B 61 , 648 (2000); 
Daniel E. Sheehy, T. P. Davis, and M. Franz, Phys. Rev. 
B 70 , 054510 (2004). 

23 Shiping Feng, Zheyu Huang, and Huaisong Zhao, Physica 
C 470 , 1968 (2010); Zheyu Huang, Huaisong Zhao, and 
Shiping Feng, Phys. Rev. B 83, 144524 (2011). 

24 Shiping Feng, Phys. Rev. B 68, 184501 (2003); Shiping 
Feng, Tianxing Ma, and Huaiming Guo, Physica C 436 , 
14 (2006). 

25 See, e.g., the review, Shiping Feng, Yu Lan, Huaisong 
Zhao, Liilin Kuang, Ling Qin, and Xixiao Ma, Int. J. Mod. 
Phys. B 29, 1530009 (2015). 

26 A. Kanigel, A. Keren, L. Patlagan, K. B. Chashka, P. King, 
and A. Amato, Phys. Rev. Lett. 92 , 257007 (2004). 

27 Y.J. Uemura, P.L.' Russo, A.T. Savici, C.R. Wiebe, G.J. 
MacDougall, G.M. Luke, M. Mochizuki, Y. Yanase, M. 
Ogata, M.L. Foo, and R.J. Cava, arXiv:cond-mat/0403031 

28 G. Baskaran, Phys. Rev. Lett. 91 , 097003 (2003); Brijesh 
Kumar and B. Sriram Shastry, Phys. Rev. B 68, 104508 
(2003); Qiang-Hua Wang, Dung-Hai Lee, and Patrick A. 
Lee, Phys. Rev. B 69 , 092504 (2004). 

29 Z. Zou and P. W. Anderson, Phys. Rev. B 37 , 627(R) 
(1988); G. Kotliar and J. Liu, Phys. Rev. B 38 , 5142(R) 
(1988). 

30 Shiping Feng, Jihong Qin, and Tianxing Ma, J. Phys.: 
Condens. Matter 16 , 343 (2004); Shiping Feng, Z. B. Su, 
and L. Yu, Phys. Rev. B 49, 2368 (1994). 

31 Bin Liu, Ying Liang, Shiping Feng, and Wei-Yeu Chen, 
Commun. Theor. Phys. 43, 1127 (2005). 

32 P. W. Anderson, Mater. Res. Bull. 8, 153 (1973); P. 
Fazekas and P. W. Anderson, Philos. Mag. 30, 423 (1974). 

33 See, e.g., the review, M. E. Zhitomirsky and A. L. Cherny¬ 
shev, Rev. Mod. Phys. 85, 219 (2013). 

34 Th. Jolicoeur and J. C. Le Guillou, Phys. Rev. B 40 , 
2727(R) (1989). 

35 B. Bernu, P. Lecheminant, C. Lhuillier, and L. Pierre, 
Phys. Rev. B 50 , 10048 (1994); P. H. Y. Li, R. F. Bishop, 
and C. E. Campbell, Phys. Rev. B 91 , 014426 (2015). 

36 Luca Capriotti, Adolfo E. Trumper, and Sandro Sorella, 
Phys. Rev. Lett. 82 , 3899 (1999). 

37 T. K. Lee and Shiping Feng, Phys. Rev. B 38 , 11809 
(1988). 

38 G. Khaliullin and P. Horsch, Phys. Rev. B 47 463 (1993). 

39 See, e.g., S. V. Tyablikov, Method in the Quantum Theory 
of Magnetism (Plenum, New York, 1967). 

40 J. Kondo and K. Yamaji, Prog. Theor. Phys. 47, 807 
(1972). 

41 N. Oeschler, R. A. Fisher, N. E. Phillips, J. E. Gordon, M.- 
L. Foo and R. J. Cava, EPL 82 , 47011 (2008); N. Oeschler, 
R. A. Fisher, N. E. Phillips, J. E. Gordon, M.-L. Foo, and 
R. J. Cava, Phys. Rev. B 78 , 054528 (2008). 



13 


42 T. K. Lee and Shiping Feng, Phys. Rev. B 41, 11110 
(1990). 

43 Carsten Honerkamp, Phys. Rev. B 68, 104510 (2003). 

44 Tsutomu Watanabe, Hisatoshi Yokoyama, Yukio Tanaka, 
Jun-ichiro Inoue, and Masao Ogata, J. Phys. Soc. Jpn. 73, 
3404 (2004). 

45 Cedric Weber, Andreas Lauchli, Frederic Mila, and Thierry 
Giamarchi, Phys. Rev. B 73, 014519 (2006). 

46 Sen Zhou and Ziqiang Wang, Phys. Rev. Lett. 100, 217002 
(2008). 

4 ' Maximilian L. Kiesel, Christian Platt, Werner Hanke, and 
Ronny Thomale, Phys. Rev. Lett. Ill, 097001 (2013). 


48 J. L. Tallon, C. Bernhard, H. Shaked, R. L. Hitterman, 
and J. D. Jorgensen, Phys. Rev. B 51, 12911(R) (1995). 

49 H. Fukuyama, H. Ebisawa, and Y. Wada, Prog. Theor. 
Phys. 42, 494 (1969); H. Fukuyama, Prog. Theor. Phys. 
42^ 1284 (1969). 

50 Mateusz Krzyzosiak, Zheyu Huang, Shiping Feng, and 
Ryszard Gonczarek, Physica C 470, 407 (2010). 

51 See, e.g., A. A. Abrikosov, Fundamentals of the Theory of 
Metals (Elsevier Science Publishers B. V., 1988). 

52 See, e.g., M. Tinkham, Introduction to Superconductivity 
(McGraw-Hill, 1996) Appendix 3. 



